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a series  of  impact  pressure  measurements  made  behind  a 
magnetically-suspended  spinning  ogive  cylinder  under 
Contract  DAAD05-74-C-0735.  This  work  was  performed  at  the 
M.I.T.  Aerophysics  Laboratory  for  the  Department  of  the 
Army,  Ballistics  Research  Laboratory,  Aberdeen  Proving  Ground, 
Md.,  under  the  technical  cognizance  of  Mr.  Anders  S.  Platou, 
Contract  Monitor.  The  bulk  of  this  material  was  part  of  a 
thesis  submitted  by  Mr.  Finn  to  the  M.I.T.  Department  of 
Aeronautics  and  Astronautics  for  the  degree  of  Master  of 
Science  in  February,  1976. 
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I. 


INTRODUCTION 


Although  the  effects  of  spin-induced  Magnus  force  on  the 
trajectories  of  axisymmetric  projectiles  were  first  ob- 
served in  the  mid-nineteenth  century,  both  theoretical  and 
experimental  research  in  the  field  continues.  The  dis- 
covery and  classical  explanation  of  the  Magnus  effect  has 
been  well  documented  in  such  reviews  as  those  of  Jacobson  (1), 
Platou  (2)  and  Swanson  (3).  Experimentation  performed  on 
large  aspect  ratio  spinning  cylinders  in  a crossflow  indicates 
a non-linear  dependance  of  the  spin-induced  force  on  spin  rate 
not  predicted  by  the  two-dimensional  theory.  Transition  has 
also  been  shown  to  grossly  affect  the  Magnus  force  in  these 
experiments  (4) . 

Extensive  theoretical  and  experimental  consideration  of  the 
Magnus  effect  has  been  applied  to  gyroscopically  stabilized 
spinning  projectiles  since  the  end  of  World  War  II.  A spinning 
body  of  revolution  at  angle  of  attack  experiences  potentially 
destabilizing  Magnus  forces  and  moments  generated  by  the  spin- 
induced  pressure  difference  on  opposing  sides  of  the  body  in 
the  crossflow  plane.  The  preponderance  of  analytical  work  has 
focused  on  the  case  of  attached  flow  under  varying  conditions 
of  spin,  angle  of  attack,  Reynolds  and  Mach  numbers.  Recent 
work  (5),  for  example,  predicts  Magnus  force  and  moment  co- 
efficients for  laminar  compressible  flow  about  several  typical 
body  geometries  by  considering  two  sources  of  the  Magnus  effect, 
these  being  circumferential  distortion  of  the  boundary  layer 
displacement  thickness  due  to  spin  and  angle  of  attack,  and  the 
variation  of  radial  pressure  gradient  due  to  centrifugal  forces 
at  the  circumference. 

Experimental  efforts  generally  encounter  several  problems 
which  are  inherent  in  the  Magnus  phenomenon.  Delicate  dynamic 
balancing  of  the  model  is  required  to  reduce  resonant  speed 
oscillations  owing  to  an  off-axis  center  of  gravity  location. 

In  addition,  forces  and  moments  measured  in  the  yaw  plane  are 
at  least  an  order  of  magnitude  less  than  the  body  normal  forces 
and  pitching  moments.  Thus,  high  sensitivity  and  repeatability 
must  be  designed  into  the  data  recording  system.  As  a result 
of  the  required  strain  gauge  sensitivity  and  the  substantial 
boundary  layer  transition  effect  on  Magnus  forces  and  moments, 
wind  tunnels  chosen  for  such  investigations  should  possess  low 
levels  of  free  stream  turbulence  and  dynamic  pressure  variation. 
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Smoke  flow  visualization  (6,7)  has  been  useful  in  delineating 
the  basic  interaction  of  the  problem  parameters  of  spin,  angle 
of  attack,  Reynold's  number  and  transition  effects.  At  small 
angles  of  attack,  low  energy  lift-induced  lee  side  vorticity 
is  shown  to  remain  embedded  in  the  boundary  layer,  thereby 
distorting  it.  For  flows  of  sufficiently  small  Reynolds 
number,  spin  induces  a skewed  laminar  boundary  layer  which 
presents  an  effectively  asymmetric  body  to  the  oncoming  flow. 

Under  laminar,  zero  spin  conditions,  higher  angles  of  attack 
cause  two  vortex  sheets  to  separate  symmetrically  from 
opposite  sides  of  the  body.  At  low  values  of  the  non- 
dimensional  spin  parameter,  PD/2V^,  the  resultant  motion  of 

the  vortex  whose  rotation  vector  is  aligned  with  that  of  the 
body  is  away  from  the  surface,  while  the  vortex  of  opposite 
rotation  remains  in  close  proximity  to  the  body  and  is  dis- 
placed in  the  direction  of  spin.  However,  flow  visualization 
reveals  that  for  high  values  of  the  spin  parameter,  the  vortex 
geometry  is  reversed.  That  is,  with  increasing  PD/2Voo,  the 

left-hand  (aligned)  vortex  is  seen  to  lose  its  identity  in 
moving  closer  to  the  body  surface  while  the  right-hand  vortex 
remains  distinct  and  is  displaced  in  the  spin  direction.  It 
has  been  suggested  in  Jacobson's  (1)  review  on  the  subject  that 
transition  from  laminar  to  turbulent  boundary  layer  conditions 
is  the  most  important  flow  condition  about  spinning  shells; 
this  is  based  on  spark  photography  (8)  results  showing  an 
asymmetric  transition  line  on  a cone  at  small  angle  of  attack 
in  a supersonic  stream. 

Recent  measurements  (9)  of  the  Magnus  force  on  a magnetically- 
suspended  spinning  five  caliber  ogive  cylinder  have  shown  that 
a sign  reversal  occurs  for  small  angles  of  attack  around  2 de- 
grees. Increasing  Reynold's  number,  addition  of  trip  devices 
and  the  introduction  of  "dummy"  sting  apparatus  into  the  base 
flow  region  all  served  to  attenuate  the  observed  reversal  effect 
and  tended  to  promote  a linear  relation  between  the  Magnus  force 
coefficient,  Cy,  and  spin  rate.  To  explain  this  it  was  postulated 

that  body  vortices  in  close  proximity  to  the  model  cause  an 
altered  pressure  distribution  which  contribute  a dominant  force 
opposite  in  direction  to  the  classical  Magnus  force  at  small 
angles  of  attack. 

In  the  present  investigation  impact  pressure  measurements  were 
made  in  the  wake  behind  a spinning  ogive  cylinder  magnetically 
suspended  in  a low  turbulence  subsonic  wind  tunnel.  Interference- 
free  support  of  the  model  in  this  facility  offers  an  excellent 
opportunity  to  study  the  wake  characteristics  corresponding  to 
the  positive  and  negative  Magnus  force  regimes. 
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To  provide  a comparative  check  of  the  drag  force  and 
corresponding  momentum  defect,  the  momentum  defect  method 
due  to  B.  M.  Jones  (10)  is  then  applied  to  calculate  profile 
drag  from  one  wake  velocity  prrf^le  behind  a non-spinning 
model.  The  drag  coefficient  so  calculated  is  then  compared 
to  force  data  from  the  magnetic  balance. 


II. 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


M.I.T.  Magnetic  Balance  Wind  Tunnel 

The  subsonic  wind  tunnel  used  in  conjunction  with  the  magnetic 
balance  system  at  the  M.I.T.  Aerophysics  Laboratory  provides 
low  turbulence  flow  and  sting-free  support  of  the  5:1  ogive 
cylinder  model.  The  open  circuit  tunnel  has  a maximum  flow 
Reynolds  number  of  five  million  per  foot,  corresponding  to 
a peak  dynamic  pressure  of  2.6  lb/in^  for  zero  blockage  con- 
ditions (11).  A free  stream  turbulence  level  of  .26%  has  been 
measured  for  Voo  = 300  ft/sec.  Dynamic  pressure  decreases 

longitudinally  in  the  test  section  by  2%  from  entrance  to  exit. 

A horizontal  buoyancy  correction  is  therefore  needed,  as  de- 
tailed in  Appendix  B.  Transverse  dynamic  pressure  variation 
does  not  exceed  ± .2%  for  distances  greater  than  or  equal  to 

one  inch  from  the  test  section  boundaries.  The  test  section  is 
a closed  octagonal  jet  made  of  plexiglass,  with  a cross-section 
area  of  32.63  in^  at  the  model  location.  The  pressure  difference 
from  two  static  taps  on  the  tunnel  contractor  is  calibrated 
versus  test  section  dynamic  pressure  providing  a means  of  setting 
flow  velocity  at  the  desired  level. 

All  pressure  measurements,  including  speed  setting,  are  made 
with  a 100  torr  full-scale  MRS  Baratron  pressure  transducer. 

The  repeatability  and  uncertainty  for  this  instrument  are  .005 
mm  Hg  and  ± .001  mm  Hg  respectively.  The  time  response  of  the 
pressure  system  was  studied  and  it  was  found  to  be  sufficiently 
short  that  it  could  be  neglected.  Model  support  in  five  degrees 
of  freedom  and  spin  about  the  horizontal  axis  is  accomplished 
with  the  magnetic  balance  system  (12) . Model  roll  is  generated 
by  a two  phase,  1200  Hz  A.C.  electromagnetic  field  which  is 
superimposed  on  the  D.C.  support  fields. 

Model  Geometry  and  Positioning 

The  wake  pressure  surveys  are  performed  in  the  flow  behind  a 
five  caliber  tangent-ogive  nose  projectile  of  one  inch  base 
diameter  (Figure  1) . The  one-pound  model  is  machined  from 
electrolytic  ingot  iron  and  is  plated  with  a copper  layer 
nominally  .003  inch  thick.  The  copper  plating  is  used  to  im- 
prove model  response  to  the  position  sensor,  decrease  system 
noise  and  improve  spin  field  coupling. 

Black  and  white  striping  was  applied  to  the  tail  end  of  the 
afterbody  to  trigger  a photocell  and  transit  focused  on  this 
portion  of  the  body.  The  amplified  output  signal  of  the 
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photocell  was  connected  to  two  digital  counters  for  monitoring 
model  spin  rate.  When  desired,  forced  boundary  layer  transition 
was  induced  by  the  application  of  #180  carborundum  grit  in  a 
one-inch  wide  strip,  beginning  one-quarter  inch  from  the  nose 
(Figure  lb) . 

In  past  work  (9)  a thin  center  body  supported  by  piano  wires 
extending  through  the  diffuser  walls  served  to  position  dummy 
stings  at  a fixed  distance  behind  a similar  five-caliber  model. 
This  stationary  body  was  used  to  support  the  impact  probe  to 
avoid  the  blocking  associated  with  a movable  pylon  mounted 
probe.  All  survey  data  in  the  wake  was  therefore  obtained  by 
translation  of  the  model  in  front  of  this  fixed  probe.  Because 
Magnus  force  acts  in  the  horizontal  plane,  most  traverses  were 
made  moving  the  model  horizontally  in  slip.  Unfortunately,  this 
model  motion  prohibits  use  of  the  helium  neon  laser  positioning 
system,  which  is  capable  of  reducing  angular  uncertainty  to 
± .02  degrees  (13). 
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Surveying  transits  were  therefore  used  to  monitor  model  trans- 
lation and  angulation  by  aligning  the  crosshairs  with  an  edge  or 
other  straight  line  feature  of  the  model.  Three  transits  were 
used  to  monitor  model  position  in  five  degrees  of  freedom 
(Figure  2) . Lift  and  drag  are  concentrated  in  one  instrument. 

A second  transit  is  aimed  at  a mirror  under  the  test  section 
which  reflects  slip  and  yaw  information.  Pitch  alone  is  viewed 
through  the  last  transit.  With  the  crosshairs  serving  as 
reference  lines,  the  model  is  brought  into  position  by  adjusting 
the  D.C.  level  fed  into  the  balance  positioning  system  from  ex- 
ternal position  command  boxes.  Model  position  uncertainty  usim 
the  transits  is  ± .001  inch  in  translation  and  + .1  degree  in 
angular  displacement.  It  is  felt  that  these  uncertainties  are 
reliable  for  all  angles  of  attack  in  the  present  investigation 
except  zero  degrees.  For  the  case  where  the  model  is  non- 
spinning and  aligned,  or  closely  aligned  with  the  flow,  unstead 
boundary  layer  effects  on  alternating  sides  of  the  model  base 
seem  to  be  the  cause  of  observed  jittery  behavior.  In  this 
case,  position  uncertainty  is  estimated  from  the  transit 
operators'  experience  as  ± .030  inch  in  translation  and  + .2  de 
grees  in  angular  displacement. 

Probe  Design 

An  impact  probe  was  designed  for  minimum  flow  disturbance  and 
minimal  flow  misalignment  error.  This  probe,  shown  in  Figure  3, 
is  constructed  of  #17  gauge  hypodermic  tubing  with  a drilled-ou* 
brass  insert  forming  the  probe  tip.  The  probe  support  assembly 
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is  shown  installed  in  the  tunnel  in  Figure  4.  The  probe  and 
the  brass  sleeve  retaining  it  comprise  a sufficiently  small 
volume  of  conductive  material  that  effects  on  magnetic  balance 
performance  are  undetectable. 

The  characteristics  of  the  probe  are  as  follows.  Reynolds 

number  based  on  sensing  hole  radius  at  the  lowest  test  velocity 

is  approximately  500;  therefore,  viscosity  effects  may  be 

neglected  (Re  > 100) . The  recovery  factor  is  100%  for  a Mach 

a 5 

number  of  .25  corresponding  to  the  high  test  ReT  = 7.50  x 10  . 

Li 

An  effective  displacement  of  the  probe  center  occurs  for 
measurements  made  in  a velocity  gradient  such  as  occurs  in 
boundary  layers  and  wakes.  This  effect  is  discussed  and  a 
correction  is  computed  in  Appendix  B. 

The  probe  is  positioned  in  the  tunnel  by  adjusting  the  piano 
wire  support  tension  (Figure  4)  while  focusing  the  various 
transits  on  a removable  probe  sheath  which  slips  over  the  brass 
sleeve  and  extends  into  the  center  of  the  test  section.  The  ex 
tension  was  necessary,  since  the  coils  of  the  magnetic  balance 
block  the  view  of  probe  tip  position.  The  probe  tip  was 
arbitrarily  located  3.0  model  base  diameters  downstream  of  the 
trailing  edge.  This  position  was  chosen  in  an  attempt  to  have 
the  probe  near  enough  to  the  model  for  good  sensitivity  but  out 
of  the  recirculation  region.  Static  pressure  measurements  in 
the  wake  of  an  axisymmetric  blunt  based  model  at  Re  = 118,000 

were  made  by  Kangovi  (14),  Figure  5.  He  found  that  wake  center 
line  static  pressure  had  returned  to  the  freestream  value  for 
X/D  = 3.0. 

At  this  axial  probe  tip  petition,  however,  the  model  could  not 
be  raised  high  enough  to  measure  wake  centerline  properties  at 
7°  angle  of  attack.  A kinked  probe  (Figure  3)  with  the  same 
tip  properties  but  with  the  tip  offset  downward  was  utilized 
for  high  incidence  testing. 

Test  Schedule 

Table  A indicates  the  test  conditions  for  the  present  investi- 
gation. Pitch  angles  of  2 and  7 degrees  were  chosen  as  a 
result  of  previous  force  measurements  (9),  which  show  these 
model  inclinations  to  produce  side  force  coefficients  of  equal 
magnitude  but  opposite  sign.  Reference  9 also  determined  that 
tests  run  at  a flow  speed  of  300  ft/sec  exhibited  significant 
non-linear  and  reversal  characteristics.  Test  velocities  of 
282  ft/sec  (Re.  = 7.50  x 105)  and  119  ft/sec  (ReT  = 3.15  x 10b) 

Li  Li 


were  chosen  to  illustrate  the  effect  of  Reynolds  number  and 
transition  on  the  Magnus  effect.  This  combination  of 
Reynolds  number  and  reduced  spin  is  shown  in  Table  I.  The 
range  of  spin  rates  was  chosen  to  utilize  the  full  capability 
of  the  magnetic  balance  for  this  particular  model  geometry. 


Table  I 


PD 

2V 


for  Run  Conditions 


Re 

Low  Spin 
High  Spin 


3.15  x 105 
. 154 
.330 


7.5  x 105 
.065 
.139 


Test  Procedure 


During  testing  the  data  was  taken  as  follows.  Ambient 
pressure,  temperature  and  relative  humidity  were  recorded 
before  and  after  each  set  of  runs.  Impact  pressure  was 
measured  as  a pressure  difference  from  atmospheric  pressure 
using  the  MKS  Baratron.  For  loss-less  flow  in  the  tunnel, 
the  pressure  difference  (p  -p  ) is  zero.  The  extent  of 

3 tl  O 

pressure  loss  was  measured  before  and  after  each  set  of  runs 

by  recording  (p  ,-p  ).  at  the  test  velocity  without  a 
1 ^at  ro  tare  1 

model.  In  addition,  static  pressure  as  detected  by  a test 

section  wall  tap  in  the  plane  of  the  probe  tip  (Figure  4) 

was  also  recorded  at  this  time. 

After  taking  the  tare  readings  the  model  was  inserted  in  the 
balance,  its  pitch  angle  and  spin  rate  were  approximately 
set  and  the  wind  was  brought  up  to  speed  by  measuring  the 
contractor  static  tap  pressure  difference.  Model  spin  during 
testing  was  held  to  within  ± 3%  of  the  specified  spin  rate. 
With  the  model  on  the  tunnel  vertical  centerline  and  its 
pitch  and  spin  set  for  the  test,  lift  position  was  adjusted 
to  obtain  the  peak  pressure  defect.  Thus,  all  traverses  in 
slip  were  obtained  on  or  near  the  vertical  wake  centerline. 
Clearly,  spin  induced  sidewash  prevents  the  exact  lift  center- 
line  position  from  being  located  using  this  method.  Since 
true  wke  characterist ics  for  any  particular  test  configura- 
tion are  not  known  a priori,  this  method  of  traversing  the 
thickest  part  of  the  wake  was  considered  the  most  practical 
and  expedient.  Additional  lift  traverses  were  also  performed 
one-half  inch  above  and  below  the  lift  centerline  position 
(Table  A) . 
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After  lift  position  was  adjusted  the  model  was  translated  in 

the  negative  slip  direction  (Table  A)  until  the  edge  of  the 

wake  was  detected,  as  indicated  by  impact  pressure,  or  until 

magnetic  support  was  degraded.  The  wake  traverse  was  then 

conducted  by  recording  (p  ,-p  ) in  1/16  inch  increments, 

at  o * 

progressing  in  the  positive  slip  direction  until  the  far  wake 
edge  was  encountered.  Checkpoints  to  determine  repeatability 
were  taken  in  the  viscous  shear  layer,  where  readings  were 
generally  stable,  in  the  wake  core  and  at  zero  tunnel  vertical 
centerline  displacement. 

In  previous  experiments  in  this  wind  tunnel,  slight  pitch  and 
yaw  flow  inclinations  have  been  detected.  Prior  to  execution 
of  the  test  schedule  an  attempt  was  made  to  remove  any  exist- 
ing flow  misalignment  relative  to  the  model  by  adjusting  pitch 
and  yaw  so  as  to  produce  an  axially  symmetric  wake.  Coarse 
traverses  were  taken  in  lift  and  slip  for  various  angles  of 
misalignment  (Figure  7)  to  predict  the  correct  model  orienta- 
tion relative  to  the  transits  for  flow  symmetry.  The  correction 
to  pitch  was  determined  to  be  - .2  degrees  and  .13  degrees  for 
yaw  (nose  rotated  counterclockwise).  Thus,  all  effective 
pitch  and  yaw  angles  (a  , 6 ) are  referenced  to  this  model 
orientation.  e 

Prior  to  taking  data  at  = 0,  2 and  7 degrees,  either  the 

straight  or  kinked  impact  probe  was  positioned  with  the  transits, 
as  discussed  in  the  section  entitled  Probe  Design.  Uncertainty 
of  the  probe  tip  in  both  slip  and  lift  was  within  + .19  inch, 
corresponding  to  an  angular  uncertainty  of  1.0  degree.  Probe 
position  remained  fixed  for  all  runs  after  the  initial  probe 
transit  positioning.  Zero  spin  profiles  in  Figures  9 and  12 
show  offsets  not  greater  than  the  above-mentioned  uncertainty 
relative  to  the  tunnel  centerline,  and  are  used  as  the  datum 
wake  position  in  deducing  subsequent  spin  and  angle  of  attack 
effects . 

The  implication  is  that  individual  wake  profiles  are  internally 
consistent.  Indeed,  one  can  compare  effects  of  changes  in 
angle  of  attack,  spin  or  Reynolds  number  on  the  shape  of  the 
profiles  but  one  cannot  compare  wake  positions  unless  the  data 
was  collected  at  the  same  tunnel  installation. 


★ 

N.B.  positive  slip  direction  places  the  probe  to  the  right 
of  the  model  when  viewed  from  the  rear. 
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■RVEYS 


The  wake  survey  ■.  >ri  discussed  in  this 

section  in  1 i . :>.♦  '■  >■  • in  addition,  the 

attenuation  of  • f et  of  the  turbulent 

boundary  lay  r addition,  the  scale 

of  pertinent  wak<  ■ established  so  that 

probes  designe  i L 1 w u i ■ work  may  be 

optimized  for  m i ■ ■ m ■ ■ . e;  , . 

The  bound, ii  ■ i model  can  be  deduced 

from  the  drag  m • t ■ ■ • (<  A Laminar  boundary 

layer  on  a clean  niod  • •xh  i ted  b-  a decreasing  drag 

coefficient  u;  to  • » 1 L number  of  approximately 

100,000.  Between  100,  oor  . d ‘ be  ol<  model  is  transi- 
tional. Above  € rippe<  mod  i appears  to  have  a 

fully-developed  turbuli  nl  boun  lai  Layei  - ' bus,  boundary 

layer  conditions  for  the  low  test  flow  velocity  are  fully 
laminar.  The  transitional  Reynolds  number  flow  regime 
represented  by  the  high  test  flow  velocity  (282  ft/sec)  is 
a favorable  test  condition  due  to  the  observed  amplification 
of  non-linear  and  reversal  effects  n past  work.  Also,  the 

boundary  layer  on  a rough . ..  • irbulent  at 

ReL  = 750,000,  th<  reby  t >rd >v  iluation  of  transition 

effects  on  wake  properties. 

It  will  be  assumed  for  the  tesl  LcU  ed  that  a 

turbulent  wake  exists.  0 rtainly  th:i.  is  ’ . io  for  instances 
when  the  boundary  layer  is  turbul,  at  tht  body'  ti  ailing 
edge.  Due  to  a velocity  profile  inflectioi  point,  however, 
the  wake  also  tend  n a laminar 

boundary  layer  persists  to  the  trailing  edge  (15). 

Total  pressure  wake  p dt  tailing  the  • fects  of  Reynolds 

number,  spin,  angle  ot  rack  arr:  transition  are  presented  in 
Figures  9 through  30.  Max  data  point  uncertainty  was 

± .1  mm  Hg , while  repeata  trying  with  test 

configuration,  was  generally  better  that  .2  mm  Hg . 

Notice  in  particular  ju : i 1 ‘ hrough  21,  which  contain  the 

offset  velocity  pr<  ■,  tor  each  particular  test  con- 
figuration, corr<  probe  lift  positions  relative 

to  the  wake  centerline.  (‘  w ' ' be  recalled  (Table  A)  that 
in  addition  to  travor  i . i<  rnw  d on  the  wake  horizontal 

centerline,  data  wa  ■ • a « rded  !m  tr  verses  one-half  inch 

above  and  below  ti  e i t • on. 
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Figure  31  represents  a typical  rear  view  of  the  flow  at  zero 
spin  for  low  (2°)  and  high  (7°)  angle  of  attack.  At  small  a 
embedded  lee  side  vortices  cause  a thickened  upper  surface 
boundary  layer  while  at  larger  values  of  a,  a symmetrical 
vortex  pair  is  shed  from  the  upper  surface. 

Evidence  of  the  separation  of  the  vortex  pair  at  7 degrees 
appears  in  the  bifurcated  centerline  profiles  of  Figure  9. 

No  such  distinct  effects  are  noticed  in  the  low  angle  of 
attack  profiles  of  Figure  12,  probably  due  to  minimal 
boundary  layer  distortion  from  the  low  intensity  embedded 
vortices . 

The  effect  of  increasing  spin  on  the  wake  geometry  for  ag  = 7° 

is  seen  in  Figures  9 through  11.  For  a =2  degrees  it  is 

shown  in  Figures  17  through  19.  This  particular  combination 
of  Reynolds  number  and  angle  of  attack,  Figure  19,  does  generate 
distinct  vortices  of  sufficient  strength  to  be  detected  by  the 
probe  located  at  X/D  = 3.0.  The  effect  of  decreased  total 
pressure  in  the  core  and  the  subsequent  displacement  with  spin 
are  clearly  observed  in  the  + 1/2"  curves  of  Figures  17  and  18. 
The  profile  amplitude  is  seen  to  be  greatly  enhanced  and  dis- 
placed at  this  higher  spin  rate,  corroborating  past  flow 
visualization  (6)  studies  and  the  boundary- layer  vortex 
geometry  hypothesis  presented  above. 

Analytical  considerations  of  the  forces  acting  on  slightly  in- 
clined spinning  bodies  of  revolution  have  as  yet  failed  to  pre- 
dict the  occurrence  of  the  reverse  Magnus  force.  The  most 
complete  theory  by  Vaughn  and  Reis  (5)  applies  to  bodies  with 
attached  laminar  boundary  layers  in  supersonic  flow  and  accounts 
for  variable  nose  shape  and  compressibility . Included  in  the 
theory  are  the  contributions  to  the  Magnus  force  of  distorted 
boundary  layer  thickness  and  radial  pressure  gradient  in  addi- 
tion to  empirical  data  on  the  effect  of  lee  side  vorticity. 

The  findings  assert  that  Magnus  forces  are  directly  proportional 
to  boundary  layer  thickness  and  local  body  radius.  Two  explicit 
non-linear  effects  are  also  delineated.  Firstly,  lee  side 
vorticity  in  the  laminar  attached  flow  is  suggested  as  the  cause 
of  side  force  non-linearity  with  angle  of  attack.  Secondly,  the 
flow  asymmetry  due  to  boundary  layer  distortion  accounts  for 
observed  Magnus  force  non-linearity  with  spin.  However,  the 
presence  of  explicit  vortex  pairs  is  not  accounted  for  by  the 
theory . 

As  angle  of  attack  is  increased  the  lee  side  vortices  gain  in 
strength,  thickening  the  upper  boundary  layer.  Sufficiently 
small  Reynolds  number  and  spin  rate  produce  a laminar  boundary 
layer  which  undergoes  separation  in  the  crossflow  plane  lower 
down  on  the  cylindrical  surface  than  for  mixed  or  fully  turbulent 
boundary  layer  conditions. 
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Profiles  of  wake  characteristics  for  a.  =7  degrees  and  P > 0 

e 

are  presented  in  Figures  20,  21  and  22.  Birtwell's  force 
measurements  for  the  same  model  (9)  have  indicated  that  at 
Re^  = 7.7  x 10-’,  this  test  configuration  results  in  the 

classical  Magnus  force.  The  lift-induced  vortex  pair  is 
clearly  seen  in  the  + 1/2"  curves  of  Figures  21  and  22.  The 
fully-developed  spin  effect  is  seen  in  the  bifurcated  pro- 
file of  Figure  21.  The  left-hand  peak  is  interpreted  as  being 
due  to  the  diminished  velocity  associated  with  the  reattached 
left-hand  vortex  while  the  separated  and  displaced  right-hand 
vortex  is  reflected  in  the  right-most  peak. 

Exactly  the  opposite  distribution  of  vortex  strength  is  ob- 
served in  the  low  Reynolds  number  test  condition  of  Figure  22, 
suggesting  a marked  change  in  the  vortex  boundary  layer  inter- 
action. 

Pressure  measurement  repeatability  was  degraded  to  a maximum 
of  1 . 5 mm  Hg  for  the  ± 1/2''  profiles  during  gritted  model  tests. 
This  is  probably  due  to  position  uncertainty  not  Ap/g  un- 
certainty since  changes  in  pressure  are  caused  both  by  model 
position  changes  and  unsteadiness  in  the  flow.  Turbulent  mix- 
ing of  the  wake  region  and  the  adjoining  laminar  free  stream 
flow  can  be  expected  to  exhibit  intermittent  turbulent  behavior. 
Thus,  unsteadiness  in  the  readings  is  attributed  to  inherent 
intermittency  in  the  viscous  shear  layer  rather  than  to  gross 
motion  of  the  wake.  Figure  27  is  a worst  case  example  of  the 
data  scatter  incurred  due  to  this  effect. 

Referring  to  Figures  28,  29  and  30  the  wake  structure  is 
essentially  preserved  with  the  onset  of  turbulent  flow. 

Measured  increases  in  Cy  due  to  turbulence,  then,  are  due  to 
enhancement  of  the  flow  asymmetry  associated  with  the  classical 
Magnus  force. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


1 


Wake  profiles  for  varying  spin  rate,  angle  of  attack,  Reynolds 
number  and  boundary  layer  conditions  were  shown  to  be  con- 
sistent with  the  flow  characteristics  relevant  to  Magnus 
testing.  < 

s 

The  impact  pressure  wake  measurements  have  demonstrated  the 
ability  to  discern  the  basic  flow  characteristics  associated 
with  an  inclined,  spinning  body  of  revolution.  Data  acquisition 
could  be  improved  with  the  insertion  of  an  oscillator  in  the 
slip  position  control  loop  to  obtain  continuous  plots  of  impact 
pressure  versus  slip  position.  Sweeping  of  the  model  in  this 
manner,  however,  involves  position  control  problems  which  could 
possibly  lead  to  extensive  redesign  of  the  existing  model 
positioning  system.  Alternatively,  a movable  probe  design 
would  allow  use  of  the  laser  system,  thereby  improving  model 
position  control  and  wake  data  repeatability.  Finally,  hot  wire 
three-dimensional  flow  angularity  wake  studies  are  recommended 
to  elucidate  the  nature  of  the  wake  beyond  the  basic  level  pre- 
sented in  this  investigation. 

Application  of  the  momentum  defect  method  (Appendix  A)  using 
a radial  distribution  of  wake  impact  pressure  resulted  in  a 
predicted  profile  drag  error  of  5.3%.  Free  stream  static 
pressure  was  assumed  constant  across  the  wake  for  this  calcula- 
tion. Subsequent  measurement  of  the  static  pressure  variation 
was  shown  to  reduce  the  difference  in  the  two  profile  drag 
calculations  to  4.8%. 
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PROFILE  DRAG  BY  MOMENTUM  DEFECT 


Theoretical  Development 

The  viscous  flow  of  a fluid  about  a solid  body  results  in  a 
net  loss  of  streamwise  momentum  equal  in  magnitude  to  the 
profile  drag  imposed  on  the  body.  Profile  drag  represents 
the  sum  of  skin  friction  drag,  owing  to  the  integral  of  body 
surface  shearing  stresses,  and  form  drag,  due  to  normal  forces 
in  the  downstream  direction.  The  methods  developed  initially 
by  Betz  and  Jones  from  consideration  of  the  momentum  equation 
provide  a means  of  calculating  the  profile  drag  from  the 
velocity  profile  of  two-dimensional  or  axially  symmetric  wake 
geometries.  A unique  feature  of  this  method  of  drag  de- 
termination is  the  essential  isolation  of  profile  and  lift 
induced  components  of  drag.  Two-dimensional  tests  (10),  for 
example,  exhibit  errors  of  the  order  of  4%  for  a lift  coefficient 
of  0.8. 

A brief  discussion  of  the  origin  and  subsequent  practical 
adaptation  of  the  momentum  defect  theory  is  presented  below. 
Consider  the  integral  form  of  the  momentum  equation  (22) 
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where  F = volume  force  per  unit  mass  of  fluid 
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The  momentum  balance  is  determined  for  steady,  viscous,  in- 
compressible fluid  flow  through  the  control  volume  defined  by 
surfaces  and  S2  in  Figure  A-l.  The  cylindrical  control 


surface  is  far  removed  from  the  body,  such  that  viscous 


forces  at  the  surface  may  be  considered  negligible.  The  net 
momentum  flux  through  S2,  comprising  the  contour  of  the  model, 

is  zero.  Clearly,  the  contribution  of  in  the  momentum 

balance  is  a drag  force  on  the  fluid  equal  and  opposite  to  that 
experienced  by  the  body.  Thus,  for  the  flow  assumptions  above, 
Equation  A. 1 applied  in  the  x direction  becomes 
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(A. 3) 


J 


(p1+pu1-p2-Pu2)dA- 


r 

T 


U'ndS  + F_ 


V 


re  F - viscous  forces  acting  at  the  control  surface  0 

subscripts  .1  and  2 denote  streamwise  flow  properties  at 
upstream  and  downstream  faces  of  respectively.  Zero 

ass  flux  through  is  expressed  as 


n dS  + 


U2~u^)  dA  = 0 


(A. 4) 


mg  the  cylindrical  control  volume  become  infinitely  long. 


u 


1 


V 

oo 


Pi  = p2 


(A. 5) 


limit,  as  the  cylindrical  faces  of  recede  from  the 

infinity,  substitution  of  the  zero  flux  condition 
1 ion  A. 4)  in  the  drag  expression  above  results  in  the 
il  relationship  between  body  profile  drag  and  wake 

ocity. 


(Voo-u2)dy  two-dimensional  wake 

(Voo-u2)rdr  axially-symmetric  wake  (A.  6) 


.■ctical  purposes  methods  utilizing  independent  initial 
ions  to  evaluate  the  above  integrals  were  developed  by 
i Jones  (23).  Both  enable  drag  prediction  from  pitot 
• wake  measurements  close  to  the  body  (within  a body 

where  static  pressure  is  substantially  different  from 
* stream  value. 

to  evaluate  the  drag  solely  over  the  velocity  profile 
on  (Figure  A-2a) , Betz  introduced  a hypothetical  flow 
it  Station  2,  which  corresponds  with  the  actual  velocity 

itution  everywhere  but  in  the  region  of  velocity  defect, 
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where  the  total  pressure  has  been  increased  to  its  free  stream 
value.  A source  located  at  the  body  is  effectively  introduced 
the  flow  as  a result  of  u!,-  The  thrust  suffered  by  V„  due 

the  source  is  incorporated  into  a momentum  balance  for  this 
hypothetical  flow.  By  considering  the  difference  between 
■ guution  A. 3 above  with  u • n = 0,  and  the  momentum  balance  for 
inviscid  source  flow,  the  drag  expression  (Equation  A. 6) 
transformed  into 


. i obviates  evaluation  of  the  integrals  beyond  the  area 
e '•  on  massed  by  the  depression  in  velocity. 

i-  method  first  applies  conditions  at  Station  3 (Figure  A-2b) 
><r  removed  from  the  body  (p^=p^)  to  Equation  A. 6. 


iying  continuity  along  a stream  tube  between  Stations  2 

r i d 


pU£dA2  = pu-^dA^ 


D 

o 


(V  -u,)dA 

00  3 


(A. 9) 


essential  assumption  of  the  method  is  that  total  pressure 
variant  between  Stations  2 and  3.  Referring  then  to 

v-.'b,  the  resultant  drag  expression,  again  non-zero  only 
■ velocity  profile  depression,  is 


) rdr 


(A. 10) 


P ~P 
^o2 


In  considering  the  repercussions  of  the  lossless  flow 
assumption,  Taylor  (16)  estimated  that  a maximum  error  of 
approximately  1.5%  would  be  incurred  in  the  drag  calculation 
for  two-dimensional  wake  measurements  performed  1/20  of  a chord 
downstream  of  an  airfoil  trailing  edge.  The  assumption  that 
1/2  p (v2  + w2)  is  negligible  at  Station  2 causes  an  additional 
inaccuracy,  resulting  in  elevated  drag  prediction.  Most  ex- 
perimental work  attempting  to  predict  drag  by  momentum  defect 
for  two-dimensional  wake  geometries  utilize  the  Jones  equation 
owing  to  the  comparative  facility  of  its  use.  No  prior  applica- 
tion of  the  momentum  defect  method  has  been  made  in  the  wake 
of  a body  of  revolution  free  of  sting  interference  effects. 

Experimental  Results  and  Analysis 

The  methods  described  in  Chapter  II  were  used  to  obtain  a radial 
distribution  of  impact  pressure  3.0  inches  downstream  of  the 
ogive  cylinder  model.  An  axially  symmetric  wake  geometry  and 
constant  free  stream  static  pressure  were  assumed.  Symmetry 
of  the  slip  traverse  data  was  obtained  by  averaging  (pat~pQ) 

about  the  profile  centerline.  Simpson's  composite  formula  for 
numerical  integration  was  used  in  evaluating  the  profile  drag 
expression  (Equation  A. 10)  in  polar  coordinate  form. 


n 

D =4rr  / p (r.)-p_(r.)  I Ip  -p  ' - fp  (r . ) -p  'I  r . Ar . (A. 

° ^ \ °2  1 2 1 °oo  00  J °2  1 | 1 1 

i=0 

As  pressure  data  was  recorded  in  1/16  inch  intervals  during 
testing,  the  radial  integration  increment,  Ar^  was  also  set  at 

1/16  inch.  Assuming  that  the  profile  drag  so  calculated  de- 
lates from  the  profile  drag  indicated  by  the  magnetic  balance, 
the  incurred  error  should  decrease  by  a factor  of  16  if  twice 
as  many  radial  intervals  are  used  (17).  A 1.6%  change  in  CD 

occurred  using  Ar^/2,  indicating  essentially  a negligible 
truncation  error. 
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The  probe  displacement  effect  described  in  Appendix  B in- 
creases the  weighting  function  of  Equation  A. 11,  thereby 
increasing  the  magnitude  of  predicted  profile  drag.  Correc- 
tions for  this  effect  as  well  as  for  solid  and  wake  blockage 
have  been  included  in  the  results  presented  here. 

As  discussed  in  the  section  entitled  Model  Geometry  and 
Positioning,  unsteady  model  motion  occurred  primarily  in  the 
lift  and  slip  planes  while  trying  to  maintain  ag  and  8g  equal 

to  zero.  Gross  asymmetries  resulted  in  the  three  impact 
pressure  profiles  recorded  for  this  model  orientation  (for 
example,  see  Figure  A-3).  Referring  to  Table  B,  the  calculated 
drag  coefficient  mean  and  standard  deviation  are  .227  and  .052 
respectively  for  these  runs.  C as  deduced  from  the  magnetic 

D0 

balance  currents  is  .196.  Considering  the  value  of  standard 
deviation,  the  resulting  agreement  in  C is  as  good  as  can  be 

0 

expected  given  the  model  position  uncertainties  at  * , 3 =0. 

G G 

Three  different  causes  of  error  were  corrected  in  an  effort  to 
improve  the  agreement  with  balance  force  measurements. 

At  zero  angle  of  attack,  high  rates  of  spin  (P  = 150  rps)  were 
found  to  have  a stabilizing  effect  on  the  model  jitters.  How- 
ever, as  would  be  expected,  a low  value  of  C (Table  B)  was 

D0 

computed  from  wake  measurements.  One  source  of  error  was  assumed 
to  be  due  to  the  effective  reduction  of  impact  pressure  defect  by 
the  spin  induced  transverse  components  of  velocity. 

Steady  non-spinning  model  behavior  for  = 2 degrees  resulted 

in  a relatively  smooth  distribution  of  impact  pressure  (Figure  A-4). 
Referring  to  Table  B,  an  error  of  5.3%  is  compatible  with  two- 
dimensional  results  for  lift  coefficients  greater  than  zero  re- 
ported in  the  past  (18).  The  validity  of  the  comparison  between 
the  momentum  defect  and  magnetic  balance  values  of  profile  drag 
coefficient  hinges  on  the  accuracy  of  the  location  of  the  wake 
vertical  centerline  during  testing.  Thus  the  error  caused  by  a 
slight  vertical  offset  in  a slip  traverse  is  incorporated  in  the 
5.3%  figure  as  well  as  the  above  effect. 


The  third  error  source  considered  was  variation  of  static 
pressure  in  the  wake.  This  static  pressure  variation  also 
sheds  an  interesting  light  on  the  appropriateness  of  the  chosen 
probe  tip  position.  Static  pressure  variation  in  the  wake  has 
been  shown  to  grossly  affect  predicted  drag  error  when  measure- 
ments are  performed  close  to  the  trailing  edge  of  an  airfoil. 
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Zwaaneveld  (19)  estimates  that  the  drag  computation  performed 
on  the  basis  of  constant  free  stream  static  pressure  (rather 
than  the  wake  centerline  or  edge  static  pressures)  generates 
the  largest  errors  when  high  wake  total  pressure  losses  are 
measured.  There  has  been  no  discussion  in  the  literature  con- 
cerning the  variation  of  the  momentum  defect  error,  incurred  by 
neglecting  static  pressure  gradient,  with  increasing  measure- 
ment distances  downstream  of  blunt  based  axisymmetric  bodies. 

Profile  drag  then  was  therefore  calculated  with  a static 
pressure  distribution  utilizing  data  measured  by  the  static 
probe,  Figure  6.  This  measured  distribution  is  shown  in 
Figure  A- 5.  Turbulence  was  not  found  to  be  of  a scale  large 
enough  to  produce  unsteady  decreases  in  static  pressure.  The 
measured  pressure  difference  (p  -p)  was  uncertain  within 

± .15%  across  the  wake.  Referring  to  Table  B,  inclusion  of 

the  static  pressure  gradient  reduced  the  error  in  C by  an  in- 

0 

crement  of  0.5%.  The  static  pressure  variation  across  the  wake 
(hen  is  seen  to  have  only  a slight  effect  on  the  profile  drag 
coefficient  predicted  by  momentum  defect.  Thus  the  chosen  probe 
tip  position  (X/D  =3.0)  is  sufficiently  far  downstream  to  avoid 
the  recirculation  zone  and  is  suggested  as  a guideline  for  the 
location  of  pressure  measuring  stations  in  similar  momentum  de- 
fect applications. 


a)  Betz  Method 


mm  Hg 


Profile  Asymmetry  due  to  Model  Position  Uncertainty 


mm  Hy 
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Lateral  Probe  Position  (inches) 


. 50x10 


. ta'ic  Pressure  Distribution  3"  Downstre.-m  of  Model  Base 


Appendix  B 


FLOW  MEASUREMENT  CORRECTIONS 


Displacement  Effect 

Impact  pressure  measurements  are  made  in  the  wake  of  an 
ogive-cylinder.  Secondary  vorticity  is  induced  in  such  a 
rotational  flow  with  the  introduction  of  a probe  causing 
an  effective  displacement  of  the  sensing  hole  center  towards 
the  region  of  higher  total  pressure.  Thus,  wake  pressures 
are  measured  at  a corrected  position  displaced  from  the  true 
geometrical  probe  centerline  in  the  direction  of  greater 
pressure.  For  right  angle  cylindrical  probes,  Reference  20, 
establishes  the  displacement  effect  correction. 


“ * °'13  + °'8  5 


AZ  = .009  inch 

where  D = probe  O.D. 

d = sensing  hole  I.D. 

AZ  = probe  centerline  displacement 


The  above  correction  is  of  particular  importance  in  the 
measurement  of  profile  drag  by  momentum  defect. 

Boundary  Interference  Effects 


A model  immersed  in  the  flow  of  an  incompressible,  viscous 
fluid  constrained  by  the  boundaries  of  a closed  octagonal 
jet  experiences  varying  interference  effects  as  a function 
of  its  position  relative  to  the  walls.  In  addition,  account 
must  be  taken  of  buoyancy  and  blockage  corrections  which  are 
independent  of  position. 

Horizontal  buoyancy  is  a continuity  effect  owing  to  boundary 
layer  growth  on  the  test  section  walls.  A longitudinal  static 
pressure  gradient  thus  necessitates  a drag  force  correction. 


FRSCKDLO  FATE  tLANK-iJOT  FIL-iiD 


In  the  present  case,  overexpansion  of  the  test  section  to 
accommodate  turbulent  boundary  layer  growth  when  in  fact 
laminar  conditions  exist  cause  a longitudinal  pressure 
increase  tending  to  draw  the  model  upstream.  The  buoyancy 
drag  correction  is 

I 


AD. 


AD. 


d£ 

d£ 


0.02  oz 


with  X = 
t = 

d£ 


body  shape  parameter 
4.8 

maximum  body  thickness 
1.0  inch 

longitudinal  static  pressure  gradient 
2.84  x 10  4 psi/inch 


Presence  of  the  model  and  its  wake  constrict  the  flow  stream- 
lines in  a closed  jet,  requiring  a correction  in  longitudinal 
flow  quantities.  A total  blockage  correction  factor  combining 
both  solid  and  wake  effects  is  estimated  as  follows. 


1 S 
4 C 

= .00602 

where  S = model  frontal  area 

C = test  section  cross  sectional  area 

and  dynamic  pressure  corrected  for  blockage  is 

% = ^00  (1  + 2 et) 
qB  = 1.012  q^ 
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When  the  model  is  experiencing  a lift  or  side  force,  the  image 
vortices  used  to  model  the  boundaries  cause  a reduction  of  in- 
duced angle  of  attack  and  induced  drag.  It  has  been  shown. 
Reference  21,  that  corrections  for  a circular  jet,  when  sub- 
stituted for  the  approximate  octagonal  jet  corrections,  admit 
a maximum  error  of  1.5%  in  the  downwash  interference  parameter, 

6.  Thus,  circular  jet  corrections  ai«_  used  here.  Referring 
to  Figure  B-l,  6 is  a function  of  model  translation  and  obeys  (22) . 


6 = 


V-(f) 


2 1 2 


The  correction  for  induced  drag  is  given  as 


AD.  = 6q 
l <»  C 


Using  typical  lift  coefficient  data  from  Reference  9 and 
assuming  a value  of  6 corresponding  to  a 1.0  inch  model  dis- 
placement, the  largest  correction  to  drag  is  .002  ounce 
(c*e  = 7 degrees)  . As  this  is  an  order  of  magnitude  less  than 

the  buoyancy  correction,  the  induced  drag  correction  is 
neglected . 

Downwash  and  sidewash  corrections  are  computed  as 


w_  = 6S 

V C L 

00 

v_  = 6S 

V C Y 


where  v,  w = sidewash  and  downwash  velocities  respectively. 

Again  basing  calculations  on  a one-inch  model  displacement, 
typical  streamline  curvature  corrections  are  shown  in  Table  C. 
It  can  be  seen  that  corrections  to  the  induced  angle  associated 
with  Magnus  force  are  an  order  of  magnitude  less  than  the  lift 
corrections  and  will  be  considered  negligible. 
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In  summary,  then,  the  following  boundary  interference  effects 
are  applied  to  the  test  data. 


qB  = 1.012  qo 


Dra„-4DB  , 

CD  - ~ q^S + S C CL 


Drag- . 02 


and  downwash  angles  are  corrected  based  on  the  value  of  6 and 
the  lift  coefficient  data  of  Reference  9. 
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MPACT  PRESSURE  SURVEY  TEST  MATRIX 
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Table  B 

PROFILE  DRAG  RESULTS 


a. 

e 

O 

P (RPS) 

Momentum  Defect 

Magnetic 

Balance 

l.1 

0 

0 

0.227 

0.196 

2. 

0 

150 

0.156 

0.196 

3. 

2 

0 

0.198 

0.209 

4. 2 

2 

0 

0.199 

0.209 

^Mean 

value  for  3 

traverses . 

static  pressure  gradient  included. 
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Table  C 


DOWNWASH  CORRECTIONS 


(°) 

C^  (Ref . 9) 

C„ (Ref . 9) 

w/V„(° 

7.0 

Jj 

0.26 

Y 

.055 

2.0 

0.11 

- 

.023 

7.0 

- 

-.013 

-.003 

2.5 

_ 

+ .013 

+ .003 
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C Profile  drag  coefficient  = Profile  drag/q^S 
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Side  force  coefficient  = Side  force/q^S 

Model  base  diameter  (in) 

Profile  drag  (lb) 

Spin  rate  (RPS) 

Non-dimensional  spin  parameter 

Static  pressure  (mm  Hg) 

Impact  pressure  (mm  Hg) 

Ambient  pressure  (mm  Hg ) 

Free  stream  dynamic  pressure  (psi) 

Radial  distance  from  wake  centerline  (in) 

Reynolds  number  based  on  probe  sensing  hole  radius 

Reynolds  number  based  on  model  length 

2 

Model  base  area  (in  ) 

Free  stream  velocity  (ft/sec) 


6.  . 

il 

0 . . 

il 


Effective  model  pitch  angle 
Effective  model  yaw  angle 
Kronecker  delta 
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